Two-dimensional turbulence of dilute polymer solutions 
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We investigate theoretically and numerically the effect of polymer additives on two-dimensional 
turbulence by means of a viscoelastic model. We provide compelling evidence that at vanishingly 
small concentrations, such that the polymers are passively transported, the probability distribution 
of polymer elongation has a power law tail: its slope is related to the statistics of finite-time Lyapunov 
exponents of the flow, in quantitative agreement with theoretical predictions. We show that at 
finite concentrations and sufficiently large elasticity the polymers react on the flow with manifold 
consequences: velocity fluctuations are drastically depleted, as observed in soap film experiments; the 
velocity statistics becomes strongly intermittent; the distribution of finite-time Lyapunov exponents 
shifts to lower values, signalling the reduction of Lagrangian chaos. 

PACS numbers: 47.27.-i 
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Since the discovery of the conspicuous drag reduction 
obtained by dissolving minute amounts of long chain 
molecules in a liquid, turbulence of dilute polymer so- 
lutions has attracted a lot of attention in view of its in- 
dustrial applications (see, e.g. Refs. [1, 2, 3]). The fluid 
mechanics of polymer solutions is appropriately described 
by viscoelastic models that are able to reproduce the rhe- 
ological behavior and many other experimental observa- 
tions [4] . For example, it has been shown by Sureshkumar 
et al. that the drag reduction effect can be captured by 
numerical simulations of the channel flow of viscoelastic 
fluids [5]. Although the parameters used in those sim- 
ulations do not match the experimental ones, the quali- 
tative agreement is remarkable, and all the hallmarks of 
the turbulent flow of polymer solutions are recovered in 
numerical experiments. 

Following this premise, it is natural to ask whether a 
two-dimensional viscoelastic model can reproduce the re- 
cent results by Amarouchenc and Kcllay [6] showing that 
the turbulent flow of soap films is spectacularly affected 
by polymer additives (see also Refs. [7, 8]). Here we show 
that this is indeed the case, and that the suppression of 
large-scale velocity fluctuations observed experimentally 
has a simple theoretical explanation. However, the influ- 
ence of polymers is not limited to the depletion of mean 
square velocity, which is a genuinely two-dimensional ef- 
fect. In the viscoelastic case we observe a strong intermit- 
tency, with exponential tails of the velocity probability 
density. As for the Lagrangian statistics, we show that 
the values of finite-time Lyapunov exponents lower signif- 
icantly upon polymer addition, which therefore reduces 
the chaoticity of the flow. These effects are expected 
to be independent of the space dimensionality, and thus 
relevant to three-dimensional turbulence as well. 

We also investigate the limit of vanishingly small poly- 
mer concentrations, in which the polymer molecules have 
no influence on the advecting flow. In this case the ve- 
locity field evolves according the two-dimensional Navier- 
Stokes equation with friction, and is therefore smooth at 



scales smaller than the injection lengthscale [9, 10]. For 
passive polymers, space dimensionality plays only a mi- 
nor role, and our system is an instance of a generic ran- 
dom smooth flow to which the theory of passive polymers 
developed by Chertkov [11] and Balkovsky et al. [12, 13] 
applies. We check this theory against our numerical re- 
sults, and find an excellent quantitative agreement. 

To describe the dynamics of a dilute polymer solution 
we adopt the linear viscoelastic model (Oldroyd-B) 

d t u+(u- V)ii = -Vp+vAu+^^-V -cr -au + f (1) 



d t a- + (u ■ V)er = (Vm) t • cr - 



■ff-(Vti)-2^-^ . (2) 



The velocity field u is incompressible, the symmetric ma- 
trix <r is the conformation tensor of polymer molecules, 
and its trace trer is a measure of their elongation [14]. 
The parameter r is the (slowest) polymer relaxation time. 
The energy source / is a large-scale random, zero-mean, 
statistically homogeneous and isotropic, solcnoidal vector 
field. The pressure term — Vp ensures incompressibility 
of the velocity field. The matrix of velocity gradients is 
defined as (Vu) !3 = diUj and 1 is the unit tensor. The 
solvent viscosity is denoted by v and r\ is the zero-shear 
contribution of polymers to the total solution viscosity 
v t = v(l + rf). The dissipative term —au models the me- 
chanical friction between the soap film and the surround- 
ing air [15], and plays a prominent role in the energy 
budget of Newtonian two-dimensional turbulence [16]. It 
should be remarked that a model that describes more 
accurately the polymer dynamics is the FENE-P model, 
which accounts for the nonlinear character of polymer 
elasticity, culminating in a finite molecular extensibil- 
ity [4]. Although here we limit ourselves to the linear 
case because it allows a simpler theoretical treatment, 
our conclusions apply to the nonlinear case as well, pro- 
vided that the maximal polymer elongation is very large 
compared to the equilibrium length. 
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Passive polymers. The effect of polymer concentration 
n is included in Eq. (1) through the parameter 77 oc n. 
In the limit rj = polymers are passively transported 
and stretched by Newtonian two-dimensional turbulence. 
The flow is driven at the largest scales and develops an 
enstrophy cascade towards the small scales, while the in- 
verse energy flux is immediately halted by friction. The 
ensuing velocity field is therefore everywhere smooth. We 
briefly recall that according to Refs.[ll, 12, 13] the statis- 
tics of stretched polymers in random smooth flows is 
expected to depend critically on the value of the Weis- 
senberg number, here defined as Wi — Xnt, where Aat 
is the maximum Lyapunov exponent of the Newtonian 
flow. At Wi < 1 the polymer molecules spend most of 
the time in a coiled state, and stretch occasionally by a 
considerable amount. The theory predicts a power law 
tail for the probability density function of trcr, i.e. the 
square polymer elongation 

p(trer) ~ (trcr)" 1 " 9 for trcr > trl . (3) 

The exponent q is related to the probability of finite-time 
Lyapunov exponents P{"), t) oc exp[— tS(j)] via the equa- 
tion L2q = 2q/r, where L,2 q — max 7 [2^7 — S(-f)] is the 
generalized Lyapunov exponent of order 2q, and 5(7) is 
the Cramer rate function (see, e.g., Ref. [17]). The con- 
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FIG. 1: Power law tail of the probability density function of 
polymer square elongation, in the passive case 7? — 0. The 
Weissenberg number is Wi = 0.4, quite below the coil-stretch 
transition. The power law (trcr)" 1 " 9 with the theoretical 
value q = 0.66 (obtained from the relation Li q = ^q/r) is 
drawn for comparison. In the inset, the corresponding Cramer 
functions 5(7). Its minimum is S(Xn) = 0, with Aat ~ 0.8. 
Data have been obtained by direct numerical simulation of 
the equation of the conformation tensor (2) by a Lagrangian 
method (see, e.g. [18] and references therein), while Eq. (1) 
has been solved in a doubly periodic box by a pseudospec- 
tral code at resolution 256 . The velocity field is driven by 
a Gaussian, homogeneous, isotropic ^-correlated in time forc- 
ing, with correlation length L w 4. The Reynolds number is 
Re — UrmsL/u w 4000. 



vexity of S("f) ensures the positivity of q for Wi < 1. 
Since the distribution of polymer elongations is not ac- 
cessible experimentally, in order to validate the theory 
it is necessary to resort to numerical simulations. Eck- 
hardt et al. in Ref. [f9] have given the first evidence 
of a power law tail for the probability distribution func- 
tion of polymer elongation in three dimensional shear 
turbulence. As shown in Fig. 1, in our two-dimensional 
simulations we observe a neat power law as well. In or- 
der to check whether the observed exponent coincides 
with the predicted one, we have also performed direct 
numerical simulations of particle trajectories, and mea- 
sured the probability distribution of finite-time Lyapunov 
exponents, thereby obtaining the expected q. The nu- 
merical result is in close agreement with theory. As the 
Weissenberg number exceeds unity, the probability distri- 
bution of the conformation tensor becomes unstationary 
and all moments ((trcr)") grow exponentially in time. 
This "coil-stretch" transition signals the breakdown of 
linear passive theory. Accounting for the nonlinear clas- 
tic modulus of polymer molecules allows to recover a sta- 
tionary statistics and to develop a consistent theory of 
passive polymers at all Weissenberg numbers [20] . In the 
following we do not pursue that approach, but we rather 
focus on a different mechanism that limits polymer elon- 
gation: the feedback of polymers on the advecting flow. 

Active polymers. When 77 > 0, polymers can affect 
significantly the velocity dynamics, provided that they 
are sufficiently elongated - a condition that is met at 
Wi > 1. This strong feedback regime is characterized 
in two dimensions by a suppression of large-scale veloc- 
ity fluctuations (see Fig. 2), an effect first observed in 
soap film experiments [6]. In Fig. 3 we present the time 
evolution of the total kinetic energy of the system, show- 
ing that after polymer injection a drastic depletion of 
kinetic energy occurs. This should be contrasted with 




FIG. 2: Snapshots of the vorticity field V x u in the New- 
tonian (left) and in the viscoelastic case with strong feedback 
(right). Notice the suppression of large-scale structures in 
the latter case. The fields are obtained by a fully dealiased 
pseudospectral simulation of eqs. (l)-(2) at resolution 256 2 . 
The viscosity is v = 1.5 • 10" 3 , r\ = 0.2, the relaxation time is 
t — 4, the energy input is F = 3.5. As customary, an artificial 
stress-diffusivity term is added to cq.(2) to prevent numeri- 
cal instabilities [21]. The corresponding Schmidt number is 
Sc = 0.25. 
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the three-dimensional case where, on the opposite, veloc- 
ity fluctuations are larger in the viscoelastic case than in 
the Newtonian one [22] . 
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FIG. 3: Dilute polymers reduce the level of velocity fluctua- 
tions J \u(x,t)\ 2 dx. Polymers are introduced in the flow at 
t = 0. In the inset, the mean square elongation j trcr(x, t) dx 
as a function of time. 

The suppression of velocity fluctuations by polymer 
additives in two-dimensional turbulence can be easily ex- 
plained in the context of the randomly driven viscoelastic 
model (l)-(2). Indeed, the average kinetic energy balance 
in the statistically stationary state reads 

F = e+^((tr<r) -trl) + cv(M 2 ) (4) 

T 

where e = ^(|Vit| 2 ) is the viscous dissipation and F is 
the average energy input, which is flow-independent for 
a Gaussian, (^-correlated random forcing /. To obtain 
eq. (4) we multiply eq. (1) by it, add to it the trace of 
cq. (2) times rjv/T, and average over space and time. 
Since in two dimensions kinetic energy flows towards 
large scales, it is mainly drained by friction, and vis- 
cous dissipation is vanishingly small in the limit of very 
large Reynolds numbers [16]. Neglecting e and observing 
that in the Newtonian case (r] = 0) the balance (4) yields 
F — a(|it| 2 )jv, we obtain 

(|«| 2 ) = (|«| 2 } J v-^((tr<r)-trl). (5) 

Since, as a consequence of incompressibility and chaotic- 
ity of the flow, it can be shown from eq. (2) that trer > 
trl, we finally have (|«| 2 ) < (|m| 2 )at, in agreement with 
numerical results. This simple energy balance argument 
can be generalized to nonlinear elastic models. As vis- 
cosity tends to zero, the average polymer elongation in- 
creases so as to compensate for the factor v in cq. (5), 
resulting in a finite effect also in the infinite Re limit. 
Since energy is essentially dissipated by linear friction, 



the depletion of (|it| 2 ) entails immediately the reduction 
of energy dissipation. The main difference between two- 
dimensional "friction reduction" and three-dimensional 
drag reduction resides in the lengthscales involved in the 
energy drain - large scales in 2D vs small scales in 3D. 

The effect of polymer additives cannot be merely rep- 
resented by a rescaling of velocity fluctuations by a given 
factor. In Fig. 4 we show the probability distribution 
of a velocity component, u x . The choice of the x direc- 
tion is immaterial by virtue of statistical isotropy. In the 
Newtonian case the distribution is remarkably close to 
the subgaussian density 7Vexp(— c|ua;| 3 ) stemming from 
the balance between forcing and nonlinear terms in the 
Navier-Stokes equation, in agreement with the predic- 
tion by Falkovich and Lebedev [23]. On the contrary, 
the distribution in the viscoelastic case is markedly su- 
pergaussian, with approximately exponential tails. This 
strong intermittency in the velocity dynamics is due to 
the alternation of quiescent low-velocity phases ruled by 
polymer feedback and bursting events where inertial non- 
linearities take over. 
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FIG. 4: Intermittency of velocity fluctuations induced by 
polymer additives. The probability density function P(u x ) 
of the velocity component u x for the Newtonian (continu- 
ous line) and for the viscoelastic case with strong feedback 
(dashed line). Same parameters as in Fig. 2. Also shown the 
distribution r(2/3)3 3/2 exp(-c|u^ | 3 )/(4ttc) with c = 2.1- 1(T 3 
(dotted line). 

Dilute polymers also alter significantly the distribu- 
tion of finite-time Lyapunov exponents -P(7, t). In Fig. 5 
the Cramer rate function S(-f) cx t^ 1 In P(j,t) is shown 
for the Newtonian and for the viscoelastic case. Since 
in the former situation the Lyapunov exponent A at is 
greater than 1/r, were the polymers passive all moments 
of elongation would grow exponentially fast. However, 
the feedback can damp stretching so effectively that after 
polymer addition A lies below 1/r. This implies a strong 
reduction of Lagrangian chaos and a decreased mixing 
efficiency. Moreover, we find that L n is smaller than n/r 
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for all n, a result which guarantees the stationarity of the 
statistics of trcr in presence of feedback, while imposing a 
less restrictive condition on £(7) than the one proposed 
in Ref. [13]. 
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FIG. 5: Finite-time Lyapunov exponents decrease in presence 
of polymers. The Cramer rate function 5(7) for the New- 
tonian (continuous line) and for the viscoelastic case with 
strong feedback (Wi = \nt = 1.6, dashed line). Viscosity 
v = 6 • 1CP 3 , relaxation time r = 2, 77 = 0.2 (dashed), n — 2 
(dotted). For sake of completeness, we also show S(-y) for a 
mild feedback case (Wi = 0.4, n — 0.2, dash-dotted line). In 
the latter case, the Lyapunov exponent is practically identical 
to the Newtonian value, and polymers affect only the right tail 
of 5(7) reducing appreciably the probability of large stretch- 
ing events 7 >• A at. 

Finally, we discuss the influence of polymer concentra- 
tion on the properties of the flow. As shown in Fig. 5, 
the Cramer functions at two very different values of r] are 
practically indistinguishable. The level of velocity fluctu- 
ations (not shown) appears to be independent of concen- 
tration as well. This property follows from the viscoelas- 
tic equations assuming that the term oc 1/r in cq. (2) 
can be neglected if polymers are substantially stretched. 
In that case the dynamics of the conformation tensor is 
invariant under rescaling by a constant factor, allowing 
to absorb r\ in the definition of cr, and making the ve- 
locity dynamics independent of concentration. This ob- 
servation poses an interesting question: can there be a 
concentration-dependent onset of friction or drag reduc- 
tion within a linear viscoelastic model ? Since polymers 
are increasingly stretched as n — > there are two alter- 
natives: either this is a singular limit, and the passive 
case is not recovered but for r\ strictly equal to zero, or 
the feedback is not uniquely ruled by polymer elongation, 
and there are other relevant mechanisms for polymer ac- 
tivity, e.g. via the creation of strong gradients of the 
conformation tensor. The limit of vanishingly small yet 
finite r\ is very demanding at the computational level and 



its investigation will require further numerical work. 
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